arXiv:1507.08935v2 [hep-ph] 17 Aug 2015 


PROCEEDINGS 

OF SCIENCE 



On the intrinsic heavy quark content of the nucieon 


Florian Lyonnet™, Aleksander Kusina^, Karol Kovari'k^,Tomas Jezo^, Fred OIness", 
logo Schienbein^ and Ji-Young 

“Southern Methodist University, Dallas, TX 75275, USA 
^Laboratoire de Physique Subatomique et de Cosmologie, 

Universite Joseph Fourier/CNRS-IN2P3/ INPG, 

53 Avenue des Martyrs, F-38026 Grenoble, France 

“Institutfur Theoretische Physik, Westfdlische Wilhelms-Universitdt Munster, 

Wilhelm-Klemm-Strafie 9, D-48149 Munster, Germany 
Universitd di Milano-Bicocca and INFN, Sezione di Milano-Bicocca, 

Piazza della Scienza 3, 20126 Milano, Italy 

E-mail: f lyonnet@mail. smu . edu, k;usina@lpsc . in2p3 . fr, 
kovarik;@particle . uni-karlsruhe . de, tomas . jezo@inib . infn . it, 
olness@mail.smu.edu, schien@lpsc.in2p3. fr, yu@lpsc.in2p3.fr 

We demonstrate that to a very good approximation the scale-evolution of the intrinsic heavy quark 
content of the nucleon is governed by non-singlet evolution equations. This allows us to analyze 
the intrinsic heavy quark distributions without having to resort to a full-fledged global analy¬ 
sis of parton distribution functions. This freedom is then exploited to model intrinsic bottom 
distributions which are so far missing in the literature in order to estimate the impact of this non- 
perturbative contribution to the bottom-quark PDF, and on parton-parton luminosities at the LHC. 
This technique can be applied to the case of intrinsic charm, albeit within the limitations outlined 
in the following. 
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1. Introduction 

Heavy quark parton distribution functions (PDFs) play an important role in several Standard 
Model (SM) and New Physics (NP) processes at the CERN Large Hadron Collider (LHC), e.g. 
see [1] for key processes involving the bottom quark PDF In the standard approach, the heavy 
quark distributions are generated radiatively, according to DGLAP evolution equations starting 
with a perturbatively calculable boundary condition. However, a purely perturbative, extrinsic, 
treatment where the heavy quarks are radiatively generated might not be adequate and there are a 
number of models that postulate a non-perturbative, intrinsic^, heavy quark component e.g. light- 
cone [2, 3] and meson cloud models [4, 5, 6]. 

Along the years, different groups have estimated the amount of intrinsic charm (IC) allowed 
in the nucleon by performing global fits, see [7] for a recent review. Among them, the two most 
recent analyses [8, 9] set significantly different limits on the allowed IC contribution highlighting 
the utility of the techniques discussed in this paper as we can freely adjust the amount of IC/IB 
contributions without having to regenerate a complete global analysis for each case. 

In this contribution, we summarize a technique which can provide IB PDFs for any generic 
non-IB PDF set, that we have introduced in [10]. Our approach exploits the fact that the intrinsic 
bottom PDF evolves (to an excellent precision) according to a standalone non-singlet evolution 
equation; that allows us to easily obtain a matched set of IB and non-IB PDFs. Note that because 
existing data entering global analyses of proton PDFs do not constrain the IB PDF, it would not be 
useful to try and obtain information on the IB content of the nucleon using a global fit. 

The rest of this paper is organized as follows. In Sec. 2, we demonstrate that to a good approx¬ 
imation the scale-evolution of the intrinsic PDF is governed by a non-singlet evolution equation, 
and construct a set of matched IC/IB PDFs. In Sec. 3, we use the IB PDFs to obtain predictions for 
parton-parton luminosities relevant at the LHC. Linally, in Sec. 4, we summarize our results and 
present conclusions. 

2. Intrinsic heavy quark PDFs 

2.1 Definition and Evolution 

In the context of a global analysis of PDLs the different parton flavors are specified via a 
boundary condifion af fhe inpuf scale /Tq which is fypically of fhe order 0’(1 GeV). Solving fhe 
DGLAP evolufion equafions wifh fhese boundary condifions allows us fo defermine fhe PDLs af 
higher scales > /Tq. A non-perfurbafive (infrinsic) heavy quark disfribufion Qi can fhen be defined 
af fhe inpuf scale /Iq as fhe difference of fhe full boundary condifion for fhe heavy quark PDL Q 
and fhe perfurbafively calculable (exfrinsic) boundary condifion Qg. 

Denofing fhe vecfor of lighf quarks as ‘q’ and fhe heavy quark disfribufion by ‘Q’ (where Q = c 
or Q = b) fhe Dokshifzer-Gribov-Lipafov-Alfarelli-Parisi (DGLAP) evolufion equafions read 


Qo + Qi = ^Qg ® g + <Si q + Pqq <Si Qo + ^QQ <Si Qi ■ 


g — Pgg ® g + Pgq 'S> q PgQ ® Qo T j 

q = Pqg®g + Pqq®q + PqQ®QQ+IfjQ^rQ{, 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


'Note that in our definition the intrinsic PDF is a twist 2 object. 
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In which we substituted Q = Qo + Qi where Qq denotes the usual radiatively generated extrinsic 
heavy quark component and Qi is the non-perturbative intrinsic heavy quark distribution.^ 

Neglecting the crossed out terms which give a tiny contribution to the evolution of the gluon 
and light quark distributions the system of evolution equations can be separated into two indepen¬ 
dent parts. For the system of gluon, light quarks and extrinsic heavy quark (g,q, Qo) one recovers 
the same evolution equations as in the standard approach without an intrinsic heavy quark compo¬ 
nent. For the intrinsic heavy quark distribution, Qi, one finds a standalone non-singlet evolution 
equation, Qi = Pqq (g) gi . 

Allowing for a small violation of the sum rule it is possible to entirely decouple the analysis of 
the intrinsic heavy quark distribution from the rest of the system. The PDFs for the gluon, the light 
quarks and the extrinsic heavy quark can be taken from a global analysis in the standard approach 
where they already saturate the momentum sum rule. On top of these PDFs the intrinsic heavy 
quark PDF can be determined in a standalone analysis using the non-singlet evolution equation. 

This induces a violation of the momentum sum rule by the term Jq dx x (Qi +Q]) which, 
however, is very small for bottom quarks.^ We will perform numerical checks of the validity of our 
approximations below. 

2.2 Modeling the boundary condition 

The BHPS model [2] predicts the following x-dependence for the intrinsic charm (IC) parton 
distribution function: 

c\ (x) = c\ (x) oc x^[6x(l -|-x)lnx-|- (1 — x)(l -|- lOx+x^)]. (2.4) 

We expect the x-shape of the intrinsic bottom distribution b\ (x) to be very similar to the one of the 
intrinsic charm distribution. Furthermore, the normalization of IB is expected to be parametrically 
suppressed with respect to IC by a factor ~ 0.1. To fix fhe freedom relafed fo fhe scale of fhe 

boundary condifion we use in fhe following fhe Same Scales boundary condifion, which remains 

2 

valid af any scale Q: b\{x,mc) = ^ci(x,njc). Finally, lef us note fhaf if would be no problem 
fo work wifh asymmefric boundary conditions, ci(x) / ci(x) and b\{x) / b\{x), as predicted for 
example by meson cloud models [5]. 

2.3 Intrinsic heavy quark PDFs from non-singlet evolution 

Eq. (2.4) is used to define the initial x-dependence at the scale of the charm mass, the nor¬ 
malization is fixed to match the one predicted by the CTEQ6.6cO fit. The IB PDE was generated 
using the Same Scales boundary conditions (see above) together with the same x-dependent input 
of Eq. (2.4). Both PDEs were then evolved according to the non-singlet evolution equation and the 
corresponding grids were produced. Note that because in our approximation, the evolution of the 
intrinsic charm and bottom PDEs is completely decoupled; the normalization of our PDEs can be 
easily changed by means of simple rescaling. 

^Strictly speaking, the decomposition of Q into Qo and Q\ is defined at the input scale where the calculable boundary 
condition for Qo is known. Consequently, Qi ■= Q — Qo is known as well. Only due to the approximations in Eqs. (2.1) 
and (2.2) it is possible to entirely decouple Qo from Qi so that the decomposition becomes meaningful at any scale. 

^It is also acceptable in case of charm provided that the allowed normalization of IC is not too big. 
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Figure 1: (a) Ratio of the CTEQ6.6 + IC 1% and CTEQ6.6cO (1%) charm distributions, (b) Ratio of the 
CTEQ6.6cO and the CTEQ6.6 gluon distributions. The results are shown as function of x for two scales, 
= 1.69 and = 10000 GeV^. 

2.4 Numerical validation 

In order to test the ideas presented above, we use the CTEQ6.6c series of intrinsic charm 
fits and in particular CTEQ6.6c0 and CTEQ6.6cl, which employ the BHPS model with 1% and 
3.5% IC probability, respectively.^ In the following we compare our approximate IC PDEs supple¬ 
mented with the central CTEQ6.6 fit, which has a radiatively generated charm distribution, with 
the CTEQ6.6cO and CTEQ6.6cl sets where IC has been obtained from global analysis. 

In Pig. 1(a), one can see that the difference between the sum cq + ci and the CTEQ6.6cO charm 
distribution is tiny at low 2^, and smaller than 5% at the higher Q^. 

The inclusion of the intrinsic charm distribution will alter the other parton distributions, most 
notably the gluon PDP. In order to gauge this effect we compare in Pig. 1(b) the gluon distribution 
from the CTEQ6.6cO analysis with the one from the standard CTEQ6.6 fit. Por small x {x <0. 1) 
the gluon PDP is not affected by the presence of a BHPS-like intrinsic charm component which 
is concentrated at large x. We note that at large-x, where most of the difference lies, the gluon 
distribution is already quite small and the uncertainty of the gluon PDP is sizable (of order of 40 - 
50% for the CTEQ6.6 set). 

We conclude that for most applications, adding a standalone intrinsic charm distribution to an 
existing standard global analysis of PDPs is internally consistent and leads only to a small error. 
Moreover, for the case of intrinsic bottom which is additionally suppressed, the accuracy of the 
approximation will be even better. Por completeness we also provided similar validation using the 
parton-parton luminosities in [10] to which we refer the reader for further details. 

3. Possible effects of IC/IB on LHC observables 

d^- 

We study the effects of both IC and IB on parton-parton luminosities, at PHC with 

^This corresponds to the values of 0.01 and 0.035 of the first moment of the charm PDF, / dxc{x), calculated at the 
input scale Qo = = 1.3 GeV. 
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Figure 2: Ratio of cc luminosities (left) and bb luminosities (right) at the LHC 14 for charm(bottom)-quark 
PDF sets with and without an intrinsic component as a function of y/x = mn / s/S. The ratio for the cc (bb) 
luminosity (solid, green line) in the left (right) hgure reaches values of 50 (17) at y/z — 0.5. In addition 
to the curves with 1% normalization (red, dashed lines) we include the results for the 3.5% normalization 
(green, solid lines) which was found to be still compatible with the current data [11]. 


s/S = 14 TeV. We use the following definition of luminosity 


dz 




1 1 0 
1 + Si j y/s Jr 


dx 


X L 


fi (x, M M) + //■ M )/<• ('^^ A: M) 


(3.1) 


where z = xiX 2 and /,• are parton distribution functions. This allows us to assess the relevance of a 
non-perturbative heavy quark component for the production of new heavy particles coupling to the 
SM fermions. 

To explore how the presence of IC and IB would affect physics observables with a non- 
negligible heavy quark initiated subprocesses, in Fig. 2 we show the ratios of luminosities for 
charm and bottom with and without an intrinsic contribution for 1% and 3.5% normalizations. Fur¬ 
thermore, since there are no experimental constraints on the IB normalization, in Fig. 2 (right), we 
also include an extreme scenario where we remove the usual ml/ml factor; thus, the first moment 
of the IB is 1% at the initial scale nic- As expected, in the case of IB the effect is smaller but for the 
bb luminosity the IB with 3.5% normalization leads to a curve which lies clearly above the error 
band of the purely perturbative result. In the extreme scenario (which is not likely but by no means 
excluded) the IB component has a big effect on the bb. 

Note that similar results for the eg and bg luminosities have been obtained but are not shown 
here, see [10]. 


4. Conclusions 

In this article, we presented a method to generate a matched IC/IB distributions for any PDF set 
without the need for a complete global re-analysis. This allows one to easily carry out a consistent 
analysis including intrinsic heavy quark effects. Because the evolution equation for the intrinsic 
heavy quarks decouples, we can freely adjust the normalization of the IC/IB PDFs. 
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For the IB, our approximation holds to a very good precision. For the IC, the error increases 
(because the IC increases), yet our method is still useful. For an IC normalization of 1-2%, the error 
is less than the PDF uncertainties at the large-x where the IC is relevant. For a larger normalization, 
although the error may be the same order as the PDF uncertainties, the IC effects also grow and 
can be separately distinguished from the case without IC. In any case, the IC/IB represents a non- 
perturbative systematic effect which should be taken into account. 

The method presented here greatly simplifies our ability to estimate the impact of the intrinsic 
heavy quark effects on the new physics searches. It can also be very useful in searching for and 
constraining the intrinsic charm and bottom components of the nucleon by itself. In particular 
in the future facilities such as an Electron Ion Collider (EIC), the Earge Hadron-Electron collider 
(EHeC), or AETER@EHC. 

The PDE sets for intrinsic charm and intrinsic bottom discussed in this analysis (1% IC, 3.5% 
IC, 1% IB, 3.5% IB) are available from the authors upon request. 
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